Bistability between different dissipative solitons in nonlinear optics 
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We report the observation of different localized structures coexisting for the same parameter 
values in an extended system. The experimental findings are carried out in a nonlinear optical 
interferometer, and are fully confirmed by numerical simulations. The existence of each kind of 
localized structure is put in relation to a corresponding delocalized pattern observed. Quantitative 
evaluation of the range of pump parameter allowing bistability between localized structures is given. 
The dependence of the phenomena on the other relevant parameters is discussed. 
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Localization of structures is a widespread phenomenon 
occurring both in conservative and in dissipative systems. 
The study of this topic aims on the one hand to an under- 
standing of the general conditions needed for the occur- 
rence of localized structures (LS); and, on the other hand, 
to explore the potentialities offered by these objects in 
view of applications, such as information transmission 
and storage. In this context, LS have been demonstrated 
in many diverse fields, including granular materials |l|, 
fluid dynamics j2| , electroconvection in liquid crystals 3 , 
chemistry and nonlinear optics . 

With specific reference to optical systems, temporal 
solitons represent a classic and very active field of re- 
search. Recently, spatial localized structures in extended 
nonlinear optical system have been widely studied as 
well 0. 

An important class of spatial LS, often referred to as 
cavity solitons, exists in many dissipative optical sys- 
tems 0. These solitons can be used as pixels for in- 
formation storage or processing 0, 0- Their existence 
has been experimentally demonstrated in several de- 
vices lint IllL I12 . i ncluding miniaturized semicon- 
ductor resonators 12, ll^ • 

In this paper we present the experimental evidence of 
bistability between different localized structures in a dis- 
sipative nonlinear optical device. The two kind of solitons 
we find differ in shape, and are separated by a discrete 
gap in their peak intensity. If used as pixels for infor- 
mation storage, these LS represent three-state variables, 
instead of the common two-state variables (" bits" ) that 
a common soliton can encode. Consequently, use of these 
new solitons would lead to an increase of log23 — 1.585 
for the information storable in a given area. 

Localized structure bistability is observed in the pres- 
ence of two modulational instabilities having different 



critical wavenumbers. The question of the interaction 
between these two instabilities has been investigated in 
[lil for semiconductor cavities. 

The existence of bistable solitons has been pre- 
dicted in the context of modified nonlinear Schroedinger 
equations0,0|, describing hght propagation in conser- 
vative systems (e.g., optical fibers in the ideal, lossless 
case). For the spatial dissipative solitons, multistabil- 
ity behaviour of either single or multi-peaked structures 
has beenpredicted numerically in nonlinear optical res- 
onators LaLauJ. To our knowledge, an experimental 
observation of bistability between different LS has not 
yet been given, neither in optics nor in other fields. 

Our experimental system is a nonlinear interferometer 
formed by a Liquid Crystal Light Valve (LCLV) with 
optical feedback. The phase ip of an initially plane laser 
beam sent onto the device evolves as Il9l 



r^ = -(V'-^o) + ^^Vi^ + /(//fa) (1) 



Ifb = la I e*'^l/2feo (^g^v' + sinc{qB ■ f)\^ (2) 



filfb) = Vsat{l - e = ) 



(3) 



where Ij, is the diffusion length of the LCLV, Ift the feed- 
back intensity on the valve, and ipo the phase retardation 
induced by the valve on the input beam in the absence 
of feedback. In (2), /q is the input intensity, B and 1 — B 
are the fractions of the field along the principal direc- 
tions of polarization of the cavity, sinc{qB • "r) represent 
the spatial filtering and * denotes a convolution product. 
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The parameters (psat, t and a represent respectively the 
saturation phase value, the response time and the sensi- 
tivity of the LCLV. Finally, I is the propagation distance 
inside the optical loop. 

Eq. (3) describes the LCLV nonlinearity. In several 
circumstances, the phenomena observed in this system 
can be understood by expanding the exponential term in 
(3) to the first order, so that a Kerr approximation is 
obtained . For the topics of interest here, it is instead 
important to take into account the saturating character 
of the nonlinearity. 

It is well known that this system displays an extreniely 
rich set of different dynamical behaviours |irl IllL 12(1 l2ll 
E3 . , and that extended as well as localized patterns 
can be formed in the transverse wavefront of the beam 
sent to the interferometer. For the present study, we set 
the system parameters at Z = 30 mm. Id ~ 18 fj,m, B ~ 
0.52, ipsat — 5, (y9o = and the spatial bandwidth qB = 
3.7. Here and in the following, the spatial frequencies 
are expressed in units of the diffractive scale qaiffr = 
2tt/ \/2\l, where A is the optical wavelength. 

In our system, one type of LS has been shown to exist 
over very broad ranges of parameters fiol Hit Esf . It has 
circular symmetry, with a bright central peak connected 
to the dark background via a series of small amplitude 
oscillations along the radial direction. Here, we report 
another completely different localized solution exists for 
the same device. The two solitons, as observed in the 
experiment, are shown in Fig. ^ Numerical simulations 
faithfully reproduce the observations. 

The most evident feature of the new LS is its triangular 
symmetry, observed both in the central peak and in the 
tails. Hence, in the following we refer to these structures 
as Triangular Solitons (TS), and to the ones with circular 
symmetry as round LS. 

A second fundamental difference between triangular 
and round LS is their peak intensity. This can be ap- 
preciated in Fig. Ic, Id. Finally the size of the central 
peak is substantially smaller for a TS than for a round 
LS, though the overall size of the two structures including 
the tails are comparable. 

The triangular soliton and the round LS shown in Fig. 
1 are observed for identical values of all the parame- 
ters, indicating bistability between the two structures. 

The existence of solitons with triangular symmetry has 
been predicted recently reported in optical svstems|26|. 
However, no bistability between RS and TS is observed 
in these systems. 

Each of these solitons can be switched on by an appro- 
priate addressing pulse. Lower intensity pulses trigger a 
round LS, higher intensity ones a TS. In these regards, 
the observed weak sensitivity to the addressing intensity 
suggests the existence of large basins of attraction for 
each soliton; hence, the noise is not expected to be a se- 
rious limit in the use of these structures for information 
storage and processing tasks. 



(a) 



(b) 




(c) 



(d) 



FIG. 1: (a), (c): Round soliton; (b), (d): Triangular soliton. 



We are now going to investigate the dependence of the 
phenomena observed on the pump intensity. The exper- 
imental state diagram of the system is shown in Fig. 2. 
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FIG. 2: Experimental state diagram of the system. Squares: 
low uniform state; empty circles: round localized structures; 
filled circles: triangular solitons. 

Starting from a very low value of input intensity and 
gradually increasing it, the lower uniform solution is the 
only state observed up to /q ~ 0.32. From this vaule 
on, round LS are observable, when addressed via proper 
initial conditions. At /q ~ 0.61 mW/cm^ the round LS 
loses its stability, and the system jumps to the triangular 
solitons branch. If, starting from a TS, the pump is de- 
creased, the structure remains stable down to Iq ~ 0.47 
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mW/cm'^, and then decays to a round LS. If, instead, 
starting from /q ~ 0.61 mW/cni^, the pump is increased, 
the TS exists up to Iq ~ 0.69 mW/cm'^, and then desta- 
biHzes via a transition to a delocahzed irregular pattern. 

The range of bistabihty between different LS is there- 
fore 0.47 < Iq < 0.61 mW/cm'^. Above ~ 2 mW/cm'^ 
the system saturates to a uniform high output intensity 
state. 

A useful insight about the origin of the LS observed can 
be gained from the observation of Fig. |21 Fig. 3(a) shows 
the homogeneous steady state (HSS) characteristic of the 
system ipQ vs alo, together with the amplitudes of the 
localized and delocalized structures found in numerical 
simulations. In Fig. 3(b) it is shown the instability region 
of the homogeneous steady state. 

In the rest of the paper we mostly discuss numerical 
results; hence it is convenient to refer to the adimensional 
parameter alo as the pump, rather than to the intensity 
Iq in physical units. In the conditions of the experiments 
here reported, a ~ 18 crn^/mW. 

The lower branch of the HSS destabilizes at alo ~ 28, 
for a critical wavenumber ~ 0.71. The critical point 
is better visible in the inset of fig. 3b, showing an en- 
larged portion of the stability balloon. The bifurcation 
from this point is strongly subcritical, and results in the 
formation of hexagons. The upper branch loses stability 
at alo ~ 20, — 1.15, with another subcritical bifur- 
cation leading to the appearence of honeycomb patterns. 

Our simulations show that both the hexagons and the 
honeycombs are unstable within most of their existence 
region. Indeed, hexagons typically decay toward a col- 
lection of spatially well separated round LS; honeycombs 
are instead destabilized either towards a collection of tri- 
angular LS, or towards a situation of developed space- 
time chaos, depending on the value of alo. In order to 
stabilize the periodic patterns and be able to track their 
amplitudes, we introduced a filter in Fourier space allow- 
ing only for the existence of the six fundamental modes 
needed to create hexagons and honeycombs, plus all their 
harmonic combination up to the cutoff frequency fixed at 
3.7. In this way we obtain the amplitudes for the delo- 
calized patterns shown in Fig. 3a. 

In the presence of the strongly subcritical patterns here 
observed, localized structures can be formed in the pa- 
rameter range where an homogeneous stable steady state 
coexists, or is not far from, a patterned state j^. These 
patterned states can be either stable or unstable for the 
parameter values at which LS exist. In the latter case, 
LS can be interpreted as residuals of some delocalized 
pattern that has lost its stability. 

In our case, the superposition of the round LS and 
of the hexagon amplitudes seen in Fig. 3a clearly indi- 
cates the relation between the localized and delocalized 
structure. As for the relation between triangular LS and 
honeycombs, the results is less obvious from the data of 
Fig. 3a. We conjecture however that such a connection 



exists also in this case. Our conjecture is based both on 
the general argument which describes the LS as local con- 
nections between a homogeneous and a patterned state; 
and, with specific reference to our case, on the triangular 
symmetry which is common to the localized structures 
and to the neighborhoods of local maxima in an honey- 
comb lattice. 

We remind that the role of an honeycomb pattern in 
gi ving rise to cavity solitons has been already pointed out 
in |l7|. for a model of semiconductor optical resonator. In 
that case, contrary to ours, dark solitons are formed by a 
dip connecting a stable bright state and a local minimum 
of the pattern. 

We concentrated our quantitative experimental inves- 
tigations on the dependence of the phenomena observed 
on the pump intensity, while the other parameters were 
kept fixed. Let us now discuss how the other parameters 
affect the scenario. An exhaustive investigation of these 
dependencies is a demanding task, which we do not face 
at the moment. We will here limit ourselves to report 
about single parameter variations, with other parameters 
kept fixed at the values used above. 

The major effect of a decrease of the system bandwidth 
qb is a narrowing of the stability range for both solitons; 
the effect is more marked for triangular solitons than for 
round LS. As a consequence, small values of qs lead to 
small ranges of soliton bistability. If is set to values 
higher than ~ 4, no substantial modifications occur with 
respect to the situation described above. 

The combined role of diffusion and diffraction 
is expressed by the adimensional parameter a = 
(2(Zd)^fco//)^/^, giving the ratio of the diffusion to diffrac- 
tion parameter. The quantity a can be viewed as an adi- 
mensional diffusion length. We are not in condition to 
have large excursion of a in the experiment. The numer- 
ics show that an increase of a with respect to the value 
used above leads to a shrinking of the stability range for 
both the round LS and the TS. The tendency is more 
marked for the triangular solitons, which have a stronger 
content of high frequency components. 

The open loop phase ipo and the weights B and 1-B of 
the feedback fields tune the S shape of the homogeneous 
steady state solution, and modify the location and critical 
wavenumbers of the modulational instabilities affecting 
the lower and upper branches. For the value of B used 
above, bistability of LS is observed for —0.4 ^ ipo ^ 0.1. 
If ipo is kept fixed at 0, the two kinds of LS are bistable 
at 0.3 < B < 0.5. 

The saturation phase (psat turns out also to be a very 
important parameter. We refer also in this case to nu- 
merical results, since the saturation parameter can be 
little varied in the experiment. If ipsat is lowered with 
respect to the value used above {(fsat = 5), the stability 
range of TS moves to higher pump intensities, and these 
structures are no more observed when if sat ^ 3.5. If fsat 
is increased, triangular solitons continue to exist as well 
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FIG. 3: Homogeneous steady state and its stability, together with localized and delocalized structure amplitude. For the values 
of parameters used, see text. In (a): Continuous line: homogeneous solution, stable parts, and hexagons bifurcating at the 
lower critical point (g — 0.71); dashed line: homogeneous solution, unstable parts; dash-dotted line: honeycombs bifurcating 
at the upper critical point (g = 1.15); empty circles: round localized structures; filled circles: triangular solitons. 



as round LS up to (psat = 18. The range of bistability has 
a broad maximum at ipsat — 6, and then shrinks, until 
disappearing at fsat — 18. Above this value, the round 
LS still exist; in correspondence of the branch formerly 
corresponding to triangular solitons, now only space-time 
chaotic delocalized patterns are observed. 

In summary, we have reported the evidence of bista- 
bility between different localized structures in a non- 
linear optical interferometer. The observation of this 
phenomenon is completely new, and certainly more ef- 
fort is required in order to obtain a full understanding 
of the sufficient and necessary conditions for its occur- 
rence. Localized structure bistability occurs within a 
complex scenario of pattern forming events; still, it is 
clear that the phenomenon is robust with respect to pa- 
rameter variations. This point, together with the gen- 
erality of the mechanism lying at its basis, suggest that 
this phenomenon could be observed in systems of differ- 
ent nature. 
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